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ABSTRACT 

We report on the identification of a near-infrared counterpart to the massive (> 
11 Mq) binary companion of pulsar J1740— 3052. An accurate celestial position of 
PSR J1740— 3052 is determined from interferometric radio observations. Adaptive op- 
tics corrected near-infrared imaging observations show a counterpart at the interfero- 
metric position of the pulsar. The counterpart has K s — 15.87±0.10 and J—K s > 0.83. 
Based on distance and absorption estimates from models of the Galactic electron and 
dust distributions these observed magnitudes are consistent with those of a main- 
sequence star as the binary companion. We argue that this counterpart is the binary 
companion to PSR J1740— 3052 and thus rule out a stellar mass black hole as the 
pulsar companion. 

Key words: binaries: general - stars: early-type - pulsars: general - pulsars: indi- 
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1 INTRODUCTION 

Among the radio pulsars that are found in binary systems, 
four belong to a subgroup where the pulsar has the charac- 
teristics of a young, isolated pulsar, but is orbiting a massive 
companion (> 3Mq) in a rather wide (50-2000 day orbital 
period) and eccentric (e > 0.5) orbit ( |Stairs|2004||van Kerk- 
|wijk et aL 2005 1 . The companion to the first of these sys- 



tems, PSRB1259-63 ( |Johnston et al.||1992[ ) was identified 
as a V = 10 early-type B2e star. PSR J0045— 7319, located 
in the SMC, is the second system and has a V — 16 B1V 
companion ( Kaspi et al.||1994 l. No information on a possi- 
ble counterpart has been published for the fourth system, 
PSRJ1638-4725 ( |Lorimer et al.|2006 |. 

The third system, PSR J1740— 3052, was discovered by 



Stairs et al.|( 2001 1 in the Parkes Multibeam Survey ( Manch- 
ester et al.||2001[ ) as a young (350 kyr) pulsar. Subsequent 
timing observations showed that the pulsar is part of a 
binary system with an orbital period of 231 d and an ec- 
centricity of 0.57. The mass function of the pulsar yielded 
a lower limit to the mass of the companion of 11 Mq. A 
candidate counterpart to PSR J1740— 3052 was identified in 
near-infrared observations. From A-band spectra, the coun- 
terpart was classified as late-type, having a spectral type 
between K5 and M3 ( Stairs et al.|200T |. This identification 
is at odds with the binary parameters which set the Roche 
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lobe radius significantly smaller than that expected for a 
late-type star at the estimated distance. Subsequent phase- 
resolved spectroscopic observations by |Tam et al.| ( |2010[ ) 
confirmed the unlikeliness of the star being the binary com- 
panion, as no significant radial velocity variations consistent 
with the expected binary orbit were found. 

The nature of the massive binary companion remains 
unclear, and could be either an early-type main-sequence 



star or a compact remnant. Though Stairs et al. (20011 



reported small variations in the dispersion measure of 
PSR J1740— 3052 near periastron, which would be consis- 
tent with a stellar companion, they could not conclusively 
rule out a stellar mass black hole as the binary companion. 

In this paper we report on interferometric radio observa- 
tions to obtain an accurate position of the PSR J1740— 3052 
and adaptive optics corrected near-infrared observations to 
search for a stellar counterpart to PSR J1740— 3052. Sec- 
tion 2 describes the interferometric observations, while Sec- 
tions details the astrometry and photometry of the near- 
infrared imaging. We discuss the results in Section 4 and 
conclude in Section 5. 



2 RADIO OBSERVATIONS 

The NRAO Very Large Array (VLA) was used to determine 
the position of the pulsar under project code AS769. The ob- 
servation was made over a 5h period starting at UTC 23:19 
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on September 3rd, 2003 using the VLA in its most extended 
A-configuration. This observation included the Pie Town an- 
tenna for even greater East- West resolution. Spectral line 
mode mode 2AC was used to provide 25 MHz of bandwidth 
centered at 1678.4 MHz on each circular polarization. The 
correlator integration time was set to 5 s. The observation 
consisted of a 600 s observation of 3C286 followed by 14 loops 
containing three sources: 90s on PMN J1751— 2524, 90s on 
PMNJ1820-2528 and 930 s on PSR J1740-3052. 

Data reduction was performed within AIPS in a stan- 
dard manner. Flux density and bandpass calibration were 
performed on 3C286. Calibrator source PMN J1751— 2524 
was used as the position reference. Its assumed posi- 
tion was taken to be aj2ooo = 17 h 51 m 51';2630, <5.j2ooo = 
-25°24'00"060 ( |Fey et al.|2004| . Phases and amplitudes de- 
termined on this source were applied to PMN J1820— 2528 
and PSR J1740— 3052. Imaging was performed using the ro- 
bust weighting factor of unity as a compromise between res- 
olution and sensitivity; the resultant beam size was 0'.'91 x 
2'.' 17 oriented 2° East of North. Wide-field imaging of the 
PSR J1740— 3052 data was performed in order to clean side- 
lobes from numerous bright sources over the entire 25' field 
of view. Most notably, the 5' diameter object known as The 



Tornado (see Gaensler et al. 2003 and references within), 
is only 8.'l Southwest of the pulsar. Due to the lack of 
short interferometer spacings, the ability to model and clean 
this source well was limited. Fortunately the pulsar lies out- 
side the worst side-lobes which emanate mostly North and 
South from the object. The pulsar was detected less than 
one beam- width from the position reported by |Stairs et al.| 
(2001) with a flux density of 0.26 ± 0.04 mjy. Despite the 
best effort to clean the side-lobes of confusing sources, the 
pulsar was found on top of residual large scale structure. 
In order to fit the position of the pulsar, linear flux density 
slices were made along the right ascension and declination 
axes through the center of the detected pulsar. Each of these 
two slices were then fit for a 1-D Gaussian and an underlying 
linear gradient. The resultant best fit position of the pulsar 
is aj2ooo = 17 h 40 m 50!001, <5 J2 ooo = -30°52'04'.'30. The la 
uncertainties on this position are 0'.'069 in right ascension 
and 0'.'21 in declination. The effect of an incorrectly mod- 
eled underlying gradient was found to result in deflections 
much smaller than the fit uncertainty so are not considered. 
The final uncertainty reported above is entirely due to mea- 
surement error. 

The 0.8 Jy VLA calibrator source PMN J1820-2528 has 
nearly the same angular separation (6? 6) from the cali- 
brator source PMN J1751— 2524 as does PSR J1740— 3052 
(6?0) and was used as an astrometric check source. 
It was imaged using the same parameters as used for 
PSR J1740— 3052. Its fit position (q J20 oo = 18 h 20 m 57:84909, 
5j20oo = -25°28'12'.'5403) differed from its cataloged posi- 
tion (QJ2000 = 18 h 20 m 57 i :8486, <5j 2 ooo = -25°28'l2'.'584;|Fey' 
|et al.|[2~004[ ) in each axis by an amount much smaller than 
the statistical uncertainty of the pulsar position. 



3 INFRARED OBSERVATIONS 

Adaptive optics corrected near-infrared imaging of the 
field containing PSR J1740— 3052 was obtained with NACO 



at the ESO Very Large Telescope on Paranal, Chile in July 
and August of 2006. A log of the observations is given in 
Table[l] Ima ges were obtained in the J, H and its-bands 
using the S27 camera, which provides a lkxlk chip with a 
0'.'027 pix~ 1 pixel scale yielding a 28" x 28" field-of-view. The 



late-type candidate counterpart of Stairs et al. ( 2001 1 was 



used as a natural guide star for the wavefront sensor which 
corrects the atmospheric distortions in the incoming wave- 
front. The N90C10 dichroic was used for K s observations, 
while the K dichroic was used for J and _ff-band imaging. 
The dithered NACO observations were used to create a sky 
frame containing the contribution of both the dark current 
and the inhomogeneities in the sky for each filter and expo- 
sure time combination. After subtraction of this sky frame, 
the science frames were registered using integer pixel offsets 
and averaged. 

Our aim was to obtain an accurate absolute astromet- 
ric calibration of the NACO images using the USNO CCD 



Astrograph Catalog (UCAC3; Zacharias et al. 2010 1 to al- 



low comparison with the interferometric pulsar position. The 
positional uncertainties of the standards in this catalog are 
20 mas for stars with R < 14 and 70 mas for R < 16, but 
with 1.3 stars per square arc-minute the star density in the 
direction of PSR J1740— 3052 is low and none of the stan- 
dards coincides with the NACO field-of-view. 

In order to allow calibration against the UCAC3 cat- 
alog, imaging observations were obtained with two further 
instruments. Narrow-band near-infrared images at 1.19 and 



2.25 fim were taken with ISAAC (Moorwood et al. 19981, 



(NAOS-CONICA; Lenzen et al. 2003 Rousset et al. 20031 



also at the VLT in June 2006. The total exposure times were 
936s and 800s respectively. This instrument has a lkxlk 
CCD with a 0"148pix" 1 pixel scale with a 2'.5 x 2.'5 field-of- 
view. The observations were corrected for dark current and 
flat-field effects and registered and averaged using standard 
routines. Finally, R, I and z-band images were obtained with 
the Wide Field Imager (WFI; |Baade et al.|1999"l ) at the 2.2 m 
ESO telescope on La Silla, Chile in August 2006. This in- 
strument is a mosaic of 8 2kx4k CCDs each with a pixel size 
of 0'.'238pix _1 and a 8' x 16' field-of-view. We only use the 
data from the CCD that contains the pulsar location. Bias- 
subtraction and flat-fielding were preformed using standard 
routines. The I and z-band images were further corrected 
for fringe frames which were constructed from the dithered 
science observations. All data reduction was done using the 
Munich Image Data Analysis System (MIDAS). 

Stars on an 8' x 8' subsection of a 45 s i?-band WFI 
image were compared to UCAC3 standards. A total of 79 
astrometric standards coincided with the image, of which 72 
were not saturated and appeared stellar and unblended. An 
astrometric solution was determined by fitting for position 
and a four parameter transformation matrix. By iteratively 
removing outliers the solution converged using 67 standards 
providing rms residuals of 0'.'073 in right ascension and 0.070 
in declination. 

This astrometric solution was transferred from the WFI 
image to the NACO images by using the 1.19pim and 
2.25 /im ISAAC images to generate tertiary and quarternary 
astrometric catalogs. Pixel positions on the ISAAC images 
were corrected for geometric distortion using the May 2005 
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Table 1. NACO Observation log. The columns give the date and 
time of the observations, the integration time t; nt , the filter used, 
the observed FWHM a of the stellar profiles and the observed 
magnitude or limit on the magnitude in this filter used. 



Date & Time (UT) 


tint (s) 


Filter 


a 


Mag. 


July f 3 


04:11-05:03 


195 X 10 


K s 


C'13 


15.87 




05:43-05:52 


50 X 5 


H 


0'.'21 


> 15.1 




05:55-06:57 


47 x 50 


J 


C'36 


> 16.7 


August 14 


00:03-00:50 


31 x 50 


K s 


(?.'25 


> 12.8 


August 21 


00:13-00:22 


50 x 5 


H 


Of/25 


> 14.1 




00:25-01:26 


240 x 10 


J 


C'39 


> 16.1 




01:29-02:00 


48 x 25 


K 3 


0'.'20 


> 13.5 



transformation provided by ESCQ The transfer of the astro- 
metric solution from the J?- band WFI image to the 1.19 fim 
ISAAC image used 56 stars, providing rms residuals of 0'.'022 
in right ascension and 0'.'023 in declination. The solution was 
then transferred to the 2.25 /im image, using 120 stars yield- 
ing rms residuals of 0'.'012 in right ascension and O'.'Oll in 
declination. Finally, 29 stars were used to transfer the solu- 
tion to the NACO K E -b&nd image obtained on July 13. The 
residuals of the transfer were 0'.'012 in right ascension and 
0'.'016 in declination. 

The total uncertainty in the astrometry of the NACO 
image is the quadratic sum of the uncertainties in the astro- 
metric solutions, which amount to 0'.'078 in right ascension 
and 0'.'076 in declination and is dominated by the uncer- 
tainties in the calibration of the WFI image by the UCAC3 
catalog. Even though multiple transfers of the astrometric 
solution are needed, it is still more accurate than the di- 
rect calibration of the 2.25 jj,m ISAAC image using 160 stars 



from the 2MASS catalog ( Skrutskie et al. 2006 1 , which yields 



residuals of O'.'ll in both coordinates. 

Combined with the uncertainty of the interferometric 
position of PSR J1740— 3052 the 95% confidence region on 
the NACO images is an ellipse with a semi-minor axis of 
0'.'26 in right ascension and a semi-major axis of 0'.'55 in 
declination. This confidence region is shown in Fig.[l] and 
is located 0'.'18 West and 0'/37 South of the late- type star, 
which itself is located just inside the 95% confidence re- 
gion. Inspection of the NACO images shows that the July 
13th, 2006 if s -band image contains a star located inside 
the error ellipse, see Fig. [2] The position of the counter- 
part in that image is a J2 ooo = 17 h 40 m 50!001 ± 0'006 and 
<5j20oo = -30 o 52'04'.'27±0'.'08, which is offset from the in- 
terferometric position of PSR J1740-3052 by 0'.'00±0'.'10 in 
right ascension and 0'.'03±0'.'22 in declination. Based on the 
stellar density of the 7v" s -band image, we estimate that the 
probability of finding a random star this close to the inter- 
ferometric position of PSR J1740— 3052 is only 0.05 per cent. 

Point-spread-function (PSF) photometry was per- 
formed on the averaged NACO images using DAOphot II 
| Stetson|p87 l running within MIDAS. An analytical Mof- 
fat profile (Moffat 19691 with a fixed exponent of ft — 2.5 



combined with an empirical look-up table was used to repre- 
sent the PSF. The parameters of the PSF were kept constant 



1 http:/ /www. eso.org/sci/facilities/paranal/ 
instruments / isaac /hist / field_distortion.html 
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Figure 1. A 8" X 8" subsection of the Kg-band images obtained 
with NACO on July 13th, 2006. The ellipse denotes the 95% con- 
fidence error region of the interferometric radio position of the 
pulsar. The bright star near the center is the late-type star which 
|Tam et al.| | |2010| showed is unlikely to be the binary companion 
to PSR J1740— 3052. This star is located on the edge of the er- 
ror region. Only a fainter star is present in the error ellipse. The 
2" X 2" box is the field shown in Fig.[2] 




Figure 2. A 2" X 2" subsection of Fig. [T] centered on the bright 
late-type star, showing the original image and with the stellar 
profile of the late-type star subtracted. Some artifacts from the 
late-type star are still present due to imperfect definition of the 
point-spread-function. The ellipse denotes the 95% confidence er- 
ror region on the interferometric position of PSR J1740— 3052. 
The counterpart to PSR J 1740— 3052 is encircled in both panels. 



over the image. Due to the low apparent stellar density on 
the NACO images, the majority of the stars on the images 
were used to determine the PSF. 

The instrumental JHK S PSF magnitudes were cali- 



brated against values from the 2MASS catalog (Skrutskie 
|et aL||2Tj56| ). Due to the small field-of-view of NACO, 4 
stars were used to determine zero-point magnitude offsets. 
No color terms were fitted for. The rms residuals of the fits 
were 0.09 mag in J-band and 0.06 mag in both H and K a . 

The counterpart to PSR J1740— 3052 is detected in the 
July 13th 7f s -band image, where it has K s = 15.87 ± 0.10. 
Limits on the magnitude of the counterpart in the other 
images were determined using simulations. For each image 
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a set of copies were generated where an artificial star of a 
certain magnitude was placed on the position of the counter- 
part. Each of these copies was then photometered using the 
same routines as the original. The 3a detection limit was set 
when the artificial star was recovered with a signal-to-noise 
of 3. The limiting magnitudes are given in TablejTJ The most 
stringent limits set J > 16.7 and H > 15.1. 



4 DISCUSSION 

We will now use the .Ks-band detection and J and ff-band 
non-detections of the counterpart and compare them with 
values predicted for the companion of PSR J1740— 3052. 

The radio timing observations by St airs et al.| ( |2001[ ) 
place several constraints on the properties of the compan- 
ion. Assuming the pulsar has a canonical mass of 1.4 Mq, 
the measurements of the orbital period and semi-major axis 
limit the mass of the companion to M c > 11 Mq. The radius 
of the companion is constrained by the radius of the Roche 
lobe at the periastron distance of the orbit. A further up- 
per limit on the radius is set by the lack of eclipses of the 
pulsar seen in the radio observations. Figure|3] shows these 
constraints. Assuming the magnetic field of the pulsar is a 
dipole, the rotational properties of the pulsar can be used to 
determine the age of the pulsar. For PSR J1740— 3052 this 
characteristic age is 350 kyr, and binary evolution dictates 
that the companion must be older than that, as the prior 



evolution of the pulsar progenitor must be involved ( Tauris 



& van den Heuvel|2006|). 



At a Galactic longitude of I = 357?81 and latitude of 
b = -0?13, PSRJ1740-3052 is located only 2?2 from the 
Galactic Center, and hence the electron and dust column 
densities are sub stantial. At a disp ersion measure (DM) of 
740.9 cm -3 pc _1 ( |Stairs et al.|200l| ) , PSR J1740-3052 is one 
of the pulsars with the largest DM. For its line-of-sight, dis- 
tance estimates from models of the Galactic electron den- 
sity distribution based on the observed DM range from 8 to 
14 kpc ([Taylor fc Cordes|1993| ) and 6.3 to llkpc ( |Cordes fc| 
Lazio|2002[ )~ 



The Cobe dust maps by Schlegel et al. ( 1998 1 set the 



maximum reddening for this line-of-sight at Eb-v = 18. Us- 
ing the Schlegel et al. (|1998|) relative extinction coefficients, 



this is equivalent to a maximum V and i^-band absorption 
of Ay < 59 and Ak < 6.7. The nearest distinct line-of- 



sight of the Galactic extinction model by [Marshall et al.| 
(2006) is within 8' of PSR J1740— 3052 and estimates K s - 
band absorption of Ak = 1.9 at a distance of 5.0 kpc, up 
to Ak = 2.5 at 10.4 kpc. This model determines the extinc- 
tion distribution by comparing 2MASS observations with 
predictions of a Galactic stellar population synthesis model 
( Robin et al.|2003|). For comparison, the Galactic extinction 



model by Drimmel et al. (2003k, where red clump stars in 



the 2MASS observations are used to derive distance depen- 
dent extinction, estimates Ak of 1.2 at 5.0 kpc, 1.7 at 10 kpc 
and 3.0 at 15 kpc. Hence, based on the observed isT s -band 
magnitude of the counterpart and a conservative estimate 
on the range of distances (d = 5 to 15 kpc) and absorptions 
(Ak = 1.2 to 3.0), the absolute A^-band magnitude of the 
counterpart falls in the range of Mk = 1.2 to —3.0. 

Stellar isochrones by |Marigo et al. | ( |2008[ l show that 
main-sequence stars have Mj — Mk = (J — K$)o > —0.25. 




20 30 
Mass (M Q ) 

Figure 3. Constraints on the properties of the binary companion 
of PSR J1740— 3052. Assuming a canonical pulsar mass of 1.4 Mq , 
the orbital parameters measured by |Stairs et al.| | |2001[ | limit the 
companion to a mass of M c > 1 1 Mq . The Roche lobe radius at 
periastron is shown with the dashed line, setting an upper limit 
to the companion radius. A further upper limit, shown with the 
dotted line, on the companion radius is set by the absence of 
eclipses of the radio signal. Stars with radii larger than this limit 
would eclipse the pulsar when it passes the line-of-sight. Stellar 
isochrones from the Padova models by |Marigo et al.| | |2008[ l are 
shown with the solid lines, ranging from zero-age main-sequence 
models to models at ages of 0.5, 1, 2, 5, 10 and 20Myr. Each 
isochrone shows the location where the star reaches absolute K B - 
band magnitudes of Mj, = —2, —3 and —4. The observed range 
of 1.2 > Mk > —3.0 shows that the companion is on the main- 
sequence. 



Hence the observed limit on J — K s > 0.83 limits the J — 
K B reddening to Ej~k > 1.08. Using the Schlegel et al. 
(|1998| extinction coefficients this sets Ak > 0.74, which 
is consistent with the estimates from Galactic absorption 
models. The observed limit on H — K B does not provide a 
better constraint. 

Figure[3] sho ws the predictions for mass and radius of 
the |Marigo et al.| ( |2008| l isochrones. They show that zero-age 
main-sequence stars with masses up to 25 Mq will have ab- 
solute ifs-band magnitudes that satisfy the observed range 
of Mk = 1-2 to —3.0. Less massive but more evolved main- 
sequence stars up to ages of 20Myr at 11 Mq also satisfy 
this range. The constraints on the mass and brightness of 
the companion of PSR J1740— 3052 show that the compan- 
ion must be on the main-sequence. Furthermore, since the 
isochrones predict that main-sequence stars more massive 
than 11 Mq have Mk < —1.5 the 7f s -band distance modu- 
lus must be larger than (m — M)k > 17.4, which suggests 
PSR J1740— 3052 is located at the far end of the distance 
ranges estimated from Galactic electron distribution mod- 
els. 
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5 CONCLUSIONS 
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